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Abstract. - The phase-ordering kinetics of the ferromagnetic 2D Ising model with uniform dis- 
order is investigated by intensive Monte Carlo simulations. Taking into account finite-time cor- 
rections to scaling, simple ageing behaviour is observed in the two-time responses and correlators. 
The dynamical exponent z and the form of the scaling functions only depend on the ratio e/T, 
where s describes the width of the distribution of the disorder. The agreement of the predictions 
of local scale-invariance generalised to z 7^ 2 for the two-time scaling functions of response and 
correlations with the numerical data provides a direct test of generalised Galilei-invariance. 
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Introduction. — A ferromagnetic system quenched 
from an initially disordered state into its coexistence phase 
with at least two equivalent equilibrium states undergoes 
phase-ordering kinetics, driven by the surface tension be- 
tween the ordered domains whose linear size grows as 
L = L(t) ~ t x l z where z is the dynamical exponent. 
For a non-conserved order-parameter it is well-known that 
z = 2, sec [1] for a review. Phase-ordering is one of the 
instances where physical ageing occurs, by which we mean 
the following properties: (i) slow (i.e. non-exponential) 
dynamics, (ii) breaking of time-translation invariancc, and 
(iii) dynamical scaling. Because of the simple algebraic 
scaling of the linear domain size L(t), two-time correla- 
tion and response functions are expected to display the 
following simple scaling forms in the ageing regime (where 
the observation time t and the waiting time s, that are 



both measured since the quench, satisfy t, s 3> f r 
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Here <f>(t, r) is the space-time-dependent order-parameter, 
whereas h(s,r) is the conjugate magnetic field (spatial 
translation-invariance will be assumed throughout this pa- 
per) and a and b are ageing exponents. The scaling func- 
tions fc,R.(y,0) ~ y~ x c,a/ z for y — > oo which defines the 
autocorrelation exponent Ac and the autoresponse expo- 
nent A/?. For phase-ordering kinetics, it is generally ad- 
mitted that 6 = and simple scaling arguments show that 
a = l/z. For an initial high-temperature state and for 
pure fcrromagncts, Ac = Xr is independent of the known 
equilibrium exponents [1-3] . More detailed information is 
contained in the form of the scaling functions fc,R.{y, u )- 
Indeed, for the phase-ordering kinetics of pure ferromag- 
nets where z = 2, it has been shown that dynamical scaling 
can be extended to a local scale-invariance (LSI) which in 
particular implies co-variance of the linear responses un- 
der transformations t 1— * tji^t + S) in time. The other 
important ingredient is the Galilei-invariance of the deter- 
ministic part of the associated stochastic Langevin equa- 
tion. Explicit predictions for the scaling functions fc,R 
follow and numerous tests have been performed in a large 
variety of models and different physical situations, see the 
recent reviews [4-6] and references therein. 

Here, we shall present tests of LSI in cases where the 
dynamical exponent z ^ 2. The foundations of the theory 
were presented some time ago in [7] and have been refor- 
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mulated recently [8] . Several tests of LSI were performed 
in systems described by linear Langevin equations in cases 
where either z = 4 [9,10] or else < z < 2 [11]. In this 
letter, we present the first test of LSI as reformulated in [8] 
with z ^ 2 in a model where the corresponding Langevin 
equation is non-linear. 

We consider a two-dimensional ferromagnetic Ising 
model with quenched disorder. The nearest-neighbour 
hamiltonian is given by [12, 13] 



"H. = — ^ ' Jij(7i(T 

(.hi) 
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The random variables are uniformly distributed over 
[1 - e/2, 1 + e/2] where < e < 2. The model has a 
second-order phase transition at a critical temperature 
T c (e) > between a paramagnetic and a (diluted) fer- 
romagnetic state. Using heat-bath dynamics such that 
the order-parameter is non-conserved and starting from a 
fully disordered initial state, phase-ordering occurs where 
the dynamical exponent is given by 



z = z(T,e) = 2 + e/T 



(3) 



This formula can be derived from phcnomcnological scal- 
ing arguments which assume that the disorder is creat- 
ing defects with logarithmically distributed barrier heights 
parametrised by the constant e [13, 14] and from field- 
theoretical studies in the Cardy-Ostlund model [15]. Sim- 
ulations of the linear domain size L(t) ~ t x l z [13] and 
of the scaling of the autorcsponse function R(t, s; 0) [16] 
also confirm eq. (3) and furthermore suggest the empirical 
identification e = e. 

Therefore, since z depends continuously on control pa- 
rameters, the disordered 2D Ising model offers a nice pos- 
sibility to test universality and especially to test LSI for 
several values of z. 

LSI-predictions. — We now state the predictions of 
LSI with an arbitrary value of z for the two-time responses 
and correlators, whose derivation will be presented else- 
where [8,17]. First, the response function reads 



R(t,s;r) = R(t,s)T^ (^^j) 
where i?(t, s) is the autoresponse function [7] 
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where a' is a further exponent and the space-time part is 
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Here a is a dimensionful, non-universal parameter and (3 
an universal exponent. Tests of eq. (5) check the time- 
dependent symmetry 1 1— > t/(jt+S). We already presented 



such tests for the 2D disordered Ising model in detail [4, 
16] . Here we shall concentrate on a detailed test of the LSI- 
prediction eq. (6) for the space-time part in the disordered 
Ising model, which is the first time that Galilei-invariance 
generalised to z ^ 2 will be checked in a non-linear model. 

In practice, it is convenient to consider the integrated re- 
sponse function (thcrmorcmancnt magnetisation), which 
is defined as 



TRM 



(t,s;r) :=h drR(t,T;r) 



(7) 



In [18], we showed that already for z = 2 MrRM(M; r ) 
may contain important finite-time corrections to scaling. 
Adapting that procedure to an arbitrary value of z, we 
find 



M TRM (t, s; r) = M cq (t -s;r)+ r A/ a 



t r 
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where ro,ri are constants and M cq (t — s;r) is the equi- 
librium contribution which vanishes for the chosen fully 
disordered initial conditions. The ageing part of the ther- 
morcmancnt magnetisation is given by 
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As for the space-time response in the pure Ising/Potts 
model quenched to T < T c [18, 19] and as for the autore- 
sponse function in the disordered Ising model [16], we find 
in many cases that sizeable corrections to scaling as de- 
scribed by the last term in eq. (8) arise and need to be 
subtracted from the simulational raw data in order to ob- 
tain a data collapse. 

Second, we consider the autocorrelation function, gener- 
alising the approach for z = 2 [6,20]. In our new approach 
to LSI, generalised Galilei-invariance with z / 2 implies 
the existence of integrals of motion involving higher pow- 
ers of the momenta. It follows that the four-point response 
function needed for the calculation of C(t, s) factorises into 
terms ^ a ^) and also depends on an 'initial' condition [8]. 
Our final prediction for the two-time autocorrelator will 
depend on the choice for that initial condition. We shall 
consider here the following two cases. 
1. As suggested by direct and straightforward comparison 
with the lattice model, we use a fully decorrelated initial 
state, with C(0,0;r) = 5{r). We find 

fc(y) = fc(y,o) (li) 

= CiJ/ Ac/z (y - l) d / 2 - 2 -Wz+2/3/ Z ( y + ^-d/z-2/3/z 
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where we have already taken into account that 6 = 0. The 
amplitude c\ remains a fitting parameter, whereas (3 is re- 
lated to \r and Ac via = Ac — A.r. 
2. From a physical point of view, we should consider 
on what time scale the scale-invariant ageing behaviour 
really sets in. When plotting the autocorrelator C(t, s) 
over against r := t — s, the data converge rapidly towards 
a plateau and only for time differences r ~ (where 
< £ < 1 is a cross-over exponent) and for s sufficiently 
large, the scaling behaviour sets in [21]. This coincides 
with the rightmost end of the plateau [22] and also with 
the time-scale on which deviations from equilibrium are 
seen in the fluctuation-dissipation ratio [21]. We there- 
fore require an estimate for the space-dependent correla- 
tor C(s + r, s;r), with t *~ S** for s — > oo and \r\ suffi- 
ciently large. Direct simulations show that in this space- 
time regime In C ~ r 2 , sec figure Id, hence we may write 
C ~ cxp(— g(r, s)r 2 ) where the function g is fixed from 
matching it with the expected scaling behaviour. From 
the above consideration scaling behaviour should set in 
for r ~ , hence we expect that g(s^, s) ~ s 2 ' z and we 
shall therefore assume a gaussian behaviour 

C(s + T,s;r) ~ cxp(-^r 2 /s 2/z ) (12) 

where v is a free parameter. This simple ansatz naturally 
generalises the equal-time correlator for phase-ordering in 
pure systems, with z = 2 [1] and may be further justified 
by generalising the derivation given in [23] for z = 2 to 
the case at hand. We then find 

fa(v) = C2V"(V - (13) 

where p is a new parameter (related to the scaling di- 
mensions of the fields involved) and ci a normalisation 
constant. Again, we have /3 = Ac — \r- In eq. (13), the 
initial correlator (12) gives rise to the factor exp(— pk 2 ) 
(/it is a constant) in the integral. Since in the application 
to the 2D disordered Ising model the dynamical expo- 
nent z > 2, from the first exponential factor in the inte- 
gral it is clear that only the small-momentum behaviour 
of the initial correlator will appreciably contribute to the 
autocorrelator scaling function. The parameter p is found 
by requiring that for y — ► 1 the autocorrelation function 
should be both nonvanishing and finite. This leads to 

p = (2/3 + d - \ c )/z (14) 

We shall see below from our numerical results that this 
relation indeed holds true in the 2D disordered Ising model 
(within the numerical error bars). 

Results. — We now compare these predictions with 
our numerical data. For the integrated response we 
simulated systems with 300 2 spins using the standard 
heat-bath algorithm. Prepared in an uncorrelated initial 



state corresponding to infinite temperatures, the system 
is quenched to the final temperature in the presence of 
a random binary field with strength h = 0.05. In order 
to obtain good statistics we averaged over typically 50000 
different runs with different initial states and different re- 
alizations of the noise. For the autocorrelation function 
we had to simulate systems containing 600 2 spins in or- 
der to avoid the appearance of finite-size effects for the 
times accessed in the simulations. It has already been no- 
ticed before [4, 16] that in disordered magnets undergoing 
phase-ordering kinetics the autocorrelation may display 
finite-time corrections to scaling, leading to the extended 
scaling form 

C(t, S ) = f c (t/ S )-s- b 'g c (t/ S ) (15) 

with b' > 0. In fact, this extended scaling form is only 
needed for dynamical exponents z > 3. Wc illustrate this 
in figure 1 for two different cases: e = 0.5 and T = 1, 
yielding z — 2.5, in panel (a), and e = 2 and T = 1, yield- 
ing z = 4, in panels (b) and (c). Whereas for z = 2.5 
one readily observes the i/s-scaling behaviour of simple 
ageing, the existence of finite-time corrections to scaling 
is obvious for z = 4. After subtracting off this correction, 
we recover the scaling of simple ageing C(t, s) = fc(t/s), 
as shown in panel (c). We list in table 2 the values of 
the exponent b' of the correction term. It is natural that 
finite-time corrections should become increasingly impor- 
tant with increasing values of z, since the domain size 
L{t) ~ t 1 ! 7, will grow more slowly for z larger and scaling 
will set in later. 

In figures 2 and 3 we show the scaling functions for the 
ageing part of the spatial thermoremanent magnetisation 
M am (t,s;r) and the autocorrelation function C(t,s) for 
several combinations of T and e. The observation of scal- 
ing of M ag0 as a function of \r\ z / s for y = t/s fixed and 
with z given by (3) with e = e gives a further confirmation 
for this dynamical exponent, see figure 2. We point out 
that the consideration of the spatio-temporal response al- 
lows for a much more demanding test of dynamical scaling 
than is possible by merely considering the autorcsponse 
alone. Our results hence strengthen the conclusions of 
a simple power-law scaling in the disordered Ising model 
reached earlier [13,14,16]. Our finding that the deviations 
from dynamical scaling seen in the raw data of the auto- 
correlation [4, 16, 24] can be explained in terms of a con- 
ventional finite-time correction makes it clear that more 
exotic proposals such as 'superageing' as proposed in [24] 
are not required. In any case, any evidence for 'superage- 
ing' would have to argue against general theoretical argu- 
ments [25] , which use basic constraints from probability to 
assert that superageing is incompatible with scaling, be- 
fore it could be considered to be conclusive. Furthermore, 
our data suggest a stronger universality in that not only 
the dynamical exponent z = 2 + e/T but also the form of 
the entire scaling functions appear to depend merely on 
the ratio e/T but not on these two control parameters sep- 
arately. We illustrate this in figure 4, where data for both 
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Fig. 1: The autocorrelation as a function of t/s for (a) e = 0.5 and T = 1, yielding z — 2.5, and (b) e = 2 and T = 1, yielding 
z — 4. Whereas in case (a) we find the scaling behaviour of simple ageing, strong corrections to scaling behaviour are seen in 
case (b). Identifying the finite-time corrections to scaling according to (15), we obtain b' = 0.075. Subtracting off the correction 
term, the scaling behaviour of simple ageing is recovered, as shown in panel (c). In (d) we show the equal-time space-dependent 
correlator C(t,t;r) for e = 0.5 and T = 1.0 for the times t = [200, 300, 500, 1000, 2000] from bottom to top. 
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Fig. 2: Scaling behaviour of the ageing part of the integrated spatio-temporal response A^trm (t, s; r), for several waiting times 
s and the values (a) e = 2.0, T = 1.0, (b) e = 1.0, T = 1.0, (c) e = 0.5, T = 0.5 and (d) e = 0.5, T = 1.0, as a function of r z /a, 
where z is given by (3) and for two fixed values of y = t/s. The full lines are the predictions of LSI, with parameters given in 
table 1. 




Fig. 3: Scaling behaviour of the autocorrelator C(t,s), for the values (a) e = 2.0, T = 1.0, (b) e = 1.0, T = 1.0, (c) e = 0.5, 
T — 0.5 and (d) e = 0.5, T = 1.0, as a function of y = t/s and for several waiting times s. The dashed lines give the prediction 
(11) with an assumed totally disordered 'initial' correlator and the full lines give the LSI-prediction (13), with the assumed 
long-ranged 'initial' correlator (12). The parameter values used are given in table 2. 
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Fig. 4: Comparison of the scaling functions of the cases e = 
1,T — 1 (filled symbols) and e — 0.5, T = 0.5 (open symbols). 
In (a) the space-time integrated response is shown for y = 2 
(upper curve) and y = 6 (lower curve). In (b) the autocorrela- 
tor is shown. The non-universal amplitudes have been changed 
appropriately. 

the spatio-temporal response and the autocorrelator in the 
two cases T = e = 1 and T = e = 0.5 are shown, which 
according to (3) should give the same dynamical exponent 
z = 3. Remarkably, not only the ageing autocorrelation 
and autorcsponse exponents indeed agree within the nu- 
merical precision of the data, but also the scaling functions 
are compatible with each other, up to normalisation (the 
non-universal amplitudes in figure 4a was changed appro- 
priately). This is a strong indication that the universal- 
ity class of the non-equilibrium kinetics of the disordered 
Ising model should not depend on T and e separately but 
rather only on the dimensionless ratio e/T. Of course, 
this observation cannot be explained in a 'superageing' 
scenario that would be incompatible with a simple scal- 
ing law Lit) ~ t x l z and with the well-accepted physical 
idea that in phase-ordering the domain size should be the 
only relevant length scale [1]. Third, when considering 
the asymptotic fall-off of the scaling functions, we find 
slightly negative values of the exponent (3 = Ac — Ar, 
see table 1. More precise data would be needed to clar- 
ify whether these results are significantly different from 
P = 0, as one obtains in pure systems [1, 11,26]. 

After these preparations, we can now consider the form 
of the scaling functions and compare our data with the 
LSI-predictions (10) and (11,13), respectively. For the 
space-time response function, the four panels in figure 2 
clearly show an excellent agreement. Since (10) does not 
depend on the initial correlations, any assumption about 
the initial correlations only enters via the values of the 
parameters a, (3 which must the same as for the autocor- 
relation. This is the first time that Galilei-invariance gen- 
eralised to z ^ 2 [8, 11] could be directly confirmed for a 
model with an underlying non-linear Langevin equation. 
It is non-trivial that a single theory is capable to repro- 
duce the shapes of the scaling functions for values of the 



dynamical exponents which vary considerably. Together 
with our previous test [16] of LSI for the autoresponse 
function in this model, this is strong evidence that LSI 
captures the essence of the dynamical scaling behaviour 
of the linear response. 

On the other hand, tests of LSI for the autocorrelation 
functions are conceptually more difficult. First, we con- 
sider a fully disordered initial state, as suggested by the 
naive comparison with the lattice model. In figure 3 we 
compare with the prediction (11), with parameter esti- 
mates listed in table 2. Clearly, one has at best a qualita- 
tive agreement and particular for smaller values of y = t/s, 
the data deviate strongly from the curve (11). Should one 
take this as a bona fide indication that LSI could not de- 
scribe the ageing of the 2D disordered Ising model? We 
have already recalled above the general arguments [21] 
which suggest that ageing only sets in for time differences 
t = t— s ~ S** , with s 3> i m icro and < C < 1. This implies 
that the form of the scaling function should rather be re- 
lated to the 'initial' correlator C(s + r, s; r). In figure Id 
we show that the leading space-dependent behaviour of 
the equal-time correlator C(s,s;r) that strongly devi- 
ates from a totally uncorrelated correlator which means 
that the system had time enough to build up strong spa- 
tially long-range correlations. We have also checked that 
the leading r-dependence is also recovered in space-time- 
dependent correlators C(t,s;r), with y = t/s = 2 or 3. 
Therefore, the assumed initial correlator eq. (12) appears 
to be in good agreement with direct numerical data for 
the asymptotic behaviour of the 'initial' correlator. That 
form is quite reminiscent to the well-known analytic forms 
found in phase-ordering kinetics of pure systems (where 
z = 2) [1] and is one of the most simple generalisations to 
z ^ 2. We see from figure 3 that the expression (13) for 
the autocorrelator C(t, s) predicted by LSI agrees nicely 
with the data and also confirm eq. (14). This is the first 
time that a consequence of the factorisation of the four- 
point response function is confirmed in a system with a 
non-linear Langevin equation. Refining the chosen form 
(12) of the initial correlator will merely lead to small cor- 
rections to (13), in particular for y — > 1. 

Summary. — 1. Taking into account strong finite- 
time corrections to scaling, we find simple ageing be- 
haviour for both the two-time responses and correlators 
of the 2D disordered Ising model. No trace of a 'super- 
ageing' behaviour is seen. 

2. Not only the dynamical exponent z = z(T,e) as given 
by (3), but also the universal form of the scaling functions 
only depends on the ratio e/T. 

3. The agreement of the space-time response predicted 
by LSI with the numerical data constitutes the first direct 
confirmation of Galilei-invariance generalised to z ^ 2 in 
a non-linear model of ageing. 

4. The precise scaling form of the autocorrelator C(t,s) 
does depend on the initial space-time correlator at the on- 
set of the ageing regime. While the usual hypothesis of a 
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Table 1: Critical exponents and parameters of LSI for the different values of e and T for the integrated response function. 



e 


T 


z 


a = a' 




a 


a 


?"0 


T\ 


0.5 


1.0 


2.5 


0.40(3) 


0.61(1) 


0.24(2) 


-0.10(5) 


0.0064(1) 


0.0025(2) 


0.5 


0.5 


3.0 


0.33(5) 


0.51(1) 


0.20(2) 


-0.06(7) 


0.0049(1) 


0.005(2) 


1.0 


1.0 


3.0 


0.33(5) 


0.51(1) 


0.20(1) 


-0.06(7) 


0.00575(2) 


-0.02(1) 


2.0 


1.0 


4.0 


0.25(2) 


0.33(1) 


0.15(2) 


-0.04(10) 


0.0365(1) 


-0.035(2) 



Table 2: Critical exponents and parameters of LSI for the different values of e and T for the autocorrelation function. 



e T 


z Xc/ Z 


Cl 


p (2/3 + rf- X c )/z v c 2 


b' 


0.5 1.0 
0.5 0.5 
1.0 1.0 
2.0 1.0 


2.5 0.570(5) 
3.0 0.490(5) 
3.0 0.490(5) 
4.0 0.320(5) 


1.51(2) 
1.48(1) 
1.48(1) 
1.62(2) 


0.14(2) 0.16(1) 0.31(2) 1.72(1) 
0.16(3) 0.14(1) 0.35(2) 1.30(1) 
0.16(3) 0.14(1) 0.36(2) 1.30(1) 
0.15(2) 0.15(1) 0.48(2) 1.12(1) 


0.075 



fully uncorrclated order-parameter does not describe the 
data, the asymptotic form eq. (12) leads to a good agree- 
ment between LSI and the data. In this way, we have con- 
firmed for the first time the factorisation of the 2n-point 
response functions, characteristic of LSI with z ^ 2 [8], in 
a non-linear model of ageing. 

The numerical results presented in this letter are strong 
evidence in favour of the existence of a local scale- 
invariance of the 'deterministic' part of the stochastic 
Langevin equation which describes the phase-ordering ki- 
netics of the 2D disordered Ising model. This is the first 
time that LSI could be confirmed for both correlators and 
responses in a non-linear model with a dynamical expo- 
nent z ^ 2. Further tests of LSI in different systems with 
z ^ 2 would be welcome. 
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